#Note: This code is for 

#platform dependent part is here. 
# If you're using Linux or OS/X, you must make 
# Replace the dev.new() call with a call to X11, 
# with X11.new() or X11(), I can't remember which. 
dev.new()   
#End of Platform dependent part.


#make 4 windows
dev.1 <- as.integer(dev.cur())
dev.new()
dev.2 <- as.integer(dev.cur())
dev.new()
dev.3 <- as.integer(dev.cur())
dev.new()
dev.4 <- as.integer(dev.cur())

#switch windows using dev.set(), passing in the desired window.
dev.set(dev.1)

mult_time <- read.csv("Royer rt_speed data.csv",header=T)
attach(mult_time)
m58 <- subset(mult_time, gender=="Male" & grade > 4)
qqnorm(m58$M_RT,main = "Royer Multiplication Response Times")
qqline(m58$M_RT)


#switch windows again to plot the histogram
dev.set(dev.2)
hist(m58$M_RT)

#compare against a true normal distribution
dev.set(dev.3)
data <-rnorm(28)
qqnorm(data)
qqline(data)

dev.set(dev.4)
qqnorm(m58$M_Speed,main = "Royer Multiplication Speeds")
qqline(m58$M_Speed)